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ABSTRACT 

Comparison  of  Predicted  Ground-Level  Airborne  Radionuclide  Concentrations 
to  Measured  Values  Resulting  from  Operation  of  the  Los  Alamos  Meson 

Physics  Facility.  (May  1993) 

William  Vandergrift  Hoak,  B.E.E.,  Villanova  University 
Chair  of  Advisory  Committee;  Dr.  Milton  E.  McLain 

A  comparison  study  of  measured  and  predicted  downwind  radionuclide 
concentrations  from  the  Los  Alamos  Meson  Physics  Facility  (LAMPF)  was 
performed.  The  radionuclide  emissions  consist  primarily  of  the  radioisotopes 
^'C,  ‘^N,  and  ‘^O.  The  gases,  vented  to  the  outside  environment  by  a  stack 
located  at  the  facility,  potentially  increase  the  radiation  exposure  at  the  facility 
boundary. 

Emission  rate,  meteorological,  and  radiation  monitoring  station  data  were 
collected  between  September  26, 1992  and  October  3, 1992.  The 
meteorological  and  emission  data  were  input  to  the  Clean  Air  Act  Assessment 
Package-1988  (CAP88-PC)  computer  code.  The  downwind  radionuclide  air 
concentrations  predicted  by  the  code  were  compared  to  the  air  concentrations 
measured  by  the  monitoring  stations.  The  code  was  found  to  slightly 
overpredict  downwind  concentrations  during  unstable  atmospheric  conditions. 
For  stable  atmospheric  conditions,  the  code  was  not  useful  for  predicting 
downwind  air  concentrations.  This  is  thought  to  be  due  to  an  underestimation  of 
hoizontal  dispersion. 
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INTRODUCTION 

The  objective  of  this  thesis  is  to  compare  theoretical  and  measured 
radionuclide  concentrations  in  air  resulting  from  the  radioactive  air  emissions 
from  the  Los  Alamos  Meson  Physics  Facility  (LAMPF).  LAMPF  is  a  half-mile 
long  linear  accelerator  that  provides  beams  of  protons  with  energies  up  to  800 
megaelectron  volts  (MeV)  at  an  average  current  of  1  milliampere.  The  proton 
beams  are  used  for  research  in  atomic  physics,  nuclear  chemistry,  radiobiology, 
and  condensed  matter  physics.  As  this  beam  of  protons  passes  through  air,  it 
activates  (transforms  stable  isotopes  into  radioactive  isotopes)  the  elements 
comprising  air.  These  activation  products  consist  primarily  of  the  radioisotopes 
and  ‘*0.  The  gases,  vented  to  the  outside  environment  by  a  stack 
located  at  the  facility,  potentially  increase  the  radiation  exposure  at  the  facility 
boundary  and  are  responsible  for  about  95  percent  of  the  offsite  dose  caused 
by  Los  Alamos  National  Laboratory.  There  is  considerable  interest  in  verifying 
the  accuracy  of  the  models  since  LAMPF  must  cease  operation  once  model 
projected  doses  offsite  approach  regulatory  limits. 

LAMPF's  highest  intensity  primary  beam  goes  to  the  meson  area  to  make  pions 
and  their  progeny  (muons  and  neutrinos).  Pions  are  created  when  positive-ion 
proton  beam  passes  thorough  two  targets  in  Area  A,  which  is  the  largest 
experimental  hall.  Made  of  carbon,  the  targets  are  contained  within  a  massive 
concrete  shield.  From  each  target,  two  beam  channels  deliver  particle  beams  to 
the  experimental  caves.  The  main  proton  beam  continues  through  the  meson 

This  thesis  follows  the  style  of  Health  Physics. 


2 

hall  to  the  Isotope  Production  and  Radiation  Effects  Facility  and  finally  is 
absorbed  in  the  beam  stop,  where  neutrons  and  neutrinos  are  produced. 

Figure  1  shows  the  accelerator  beam  line  and  the  exhaust  system.  The 
proton  beam  is  directed  into  two  targets  (A-1  and  A-2)  which  are  used  to 
produce  subatomic  particles  for  use  in  various  experiments.  The  remaining 
beam  is  directed  into  a  beam  stop  (Shire  1992). 

Los  Alamos  National  Laboratory  is  required  by  the  Environmental 
Protection  Agency  (EPA)  to  maintain  a  dose  of  less  than  0.1  milliSieverts  (mSv) 
at  the  boundary  of  the  laboratory  (i.e.,  the  maximally  exposed  member  of  the 
general  public).  In  October  1969  the  EPA  issued  final  rules  for  radionuclide 
emissions  to  air  under  40  CFR  Part  61 ,  National  Emission  Standards  for 
Hazardous  Air  Pollutants  (NESHAPS).  Emission  monitoring  and  compliance 
procedures  for  Department  of  Energy  (DOE)  facilities  (40  CFR  61 .93(a))  require 
the  use  of  Clean  Air  Act  Assessment  Package-1988  (CAP88-PC)  or  AIRDOS- 
PC  computer  models,  or  other  approved  procedures,  to  calculate  effective  dose 
equivalents  to  members  of  the  public. 

The  CAP88'PC  model  uses  a  modified  Gaussian  plume  equation  to 
estimate  the  average  dispersion  of  radionuclides  released  from  up  to  six 
sources.  The  sources  may  be  either  elevated  stacks,  such  as  a  smokestack,  or 
uniform  area  sources.  Assessments  are  done  for  a  circular  grid  of  distances 
and  directions  of  30  kilometers  around  the  facility.  Because  of  the  simplicity  of 
the  model  and  the  conservative  assumptions  used  in  CAP88-PC,  the  EPA 
claims  that  the  doses  estimations  are  good  within  a  factor  of  two,  normally 
overprojecting  doses  for  most  situtations  (EPA  1992). 


LOS  ALAMOS  CANYON 


RADIOACTIVE 
AIR  EXHAUST 
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Because  of  the  above  short  comings  of  the  model,  doses  may  be  significantly 
overestimated  by  CAP88-PC  model.  Additionally,  it  assumes  a  that  the 
topology  of  the  area  in  question  is  a  flat  plane.  LAMPF  is  located  on  the 
Pajarito  Plateau  on  the  eastern  flanks  of  the  Jemez  Mountains  in  Northern  New 
Mexico.  There  are  numerous  alternating  finger  mesas  and  canyons  running 
along  the  slope  line  of  the  plateau.  The  canyons  are  50-1 00  m  deep  and  1 00- 
200  m  wide.  LAMPF  is  located  on  a  mesa  top  just  south  of  Los  Alamos  Canyon. 
The  town  of  Los  Alamos  is  located  on  the  north  side  of  the  mesa.  This  unique 
topography  lends  itself  to  meteorological  conditions  which  are  not  well 
addressed  by  the  model.  The  accelerator  in  relation  to  its  surrounding  is  shown 
in  Fig  2. 


ig  2.  LAMPF  and  radiation  monitoring  stations 
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BACKGROUND 


Bowen  (Bowen  et  al.  1987)  addressed  the  accuracy  of  standard 
atmospheric  models  using  pressurized  ion  chambers  for  population  dose 
estimates  in  1985.  Both  monitoring  and  modeling  of  external  radiation  levels  at 
three  locations  during  1985  were  studied.  Predicted  daily  levels  using  a  simple 
(infinite  cloud)  Gaussian  model  were  compared  with  daily  measured  values. 
Predicted  values  using  a  more  complex  finite  model  were  also  examined. 
Average,  hourly  predicted  and  measured  external  radiation  levels  were 
compared  with  each  other  and  in  relation  to  average  winds.  Accuracy  of  the 
model  was  tested  for  different  times  of  the  day,  wind  speeds,  horizontal  wind, 
and  vertical  wind.  Cumulative  external  radiation  during  an  86-day  period,  as 
measured  by  the  ion  chambers  was  compared  with  both  the  sum  of  24-hour 
estimates  and  the  seasonal  form  of  the  model  using  average  winds.  Annual 
model  predictions  are  also  compared  with  TLD  measurements. 

The  simple  infinite  cloud  model  was  found  to  overestimate  the  population 
dose  by  over  50%  using  the  ion  chambers  as  a  reference.  The  more  complex 
finite  cloud  model  over  estimated  the  dose  by  15%.  Bowen  found  that  the 
models  were  particularly  poor  for  predicting  daytime  radiation  levels.  No 
consideration  was  given  to  amount  of  radionuclides  that  are  mixed  down  into 
the  Los  Alamos  Canyon  in  transport,  or  to  incorporating  wind  shear  into  the 
calculation.  The  canyon  effect  on  diffusing  radionuclides  was  not  compared 
using  radioactivity  levels  resulting  from  transport  over  the  mesa  with  transport 
over  the  canyon. 

These  models  were  adequate  at  the  time,  because  the  dose  limits  were 
0.25  mSv  and  the  calculated  doses  were  well  below  this  level.  With  the  recent 
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change  in  the  dose  limit  to  0. 1  mSv,  the  current  models  could  under  worse  case 
conditions  estimate  doses  in  excess  of  the  0.1  mSv  limit,  requiring  the  shutdown 
of  LAMPF  for  the  remainder  of  the  year. 

Meteorological  Stations 

Meteorological  data  for  this  study  was  collected  from  an  instrumented 
meteorological  tower  located  at  the  accelerator  site  (technical  area  53).  Data 
on  horizontal  wind  direction  and  speed  are  taken  at  heights  of  11 .5,  23,  and  46 
m.  Vertical  velocity,  temperature,  solar  radiation,  relative  humidity,  and  rainfall 
are  also  recorded.  Data  is  continuously  logged  by  microprocessors  at  the  tower 
and  averaged  over  1 5  minute  intervals.  The  1 5  minute  data  is  retrieved  by 
computer  and  stored  on  the  Los  Alamos  central  computer  system  for  access  by 
end  users. 

Radiation  Monitoring  Stations 

To  more  accurately  measure  radionuclide  activities  and  compositions,  three 
gamma  radiation  spectrometers  have  been  setup  at  the  north  boundary  of  the 
laboratory.  These  computer  controlled  instruments  have  been  remotely 
collecting  data  on  hourly  intervals  since  the  last  week  of  August  1 992.  Data 
from  these  spectrometers  has  significant  advantages  over  the  pressurized  ion 
chambers  for  purposes  of  this  project.  The  natural  background  radiation  levels 
at  the  site  boundary  range  from  one  to  two  mSv  Without  spectral  data,  it  is 
very  difficult  to  discern  doses  on  the  order  of  the  0.01  mSv  y-’,  the  level  of 
accuracy  required  to  demonstrate  compliance  with  the  law. 

The  monitoring  stations  are  operated  by  the  Environmental  Management 
Division.  Only  two  stations  were  operational  during  the  measurement  period. 
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The  first  station,  referred  to  as  the  center  station,  is  located  742.8  m  from  the 
stack  along  a  22.5°  radial  from  true  north.  The  second  station,  referred  to  as 
the  west  station,  is  located  699.8  m  from  the  stack  along  a  0°  radial  (true  north). 
Each  station  is  equipped  with  a  p-type  hyperpure  germanium  (HpGe)  gamma 
spectrometer  detector  (EG&G  Ortec  model  GEM,  40%  nominal  efficiency)  with 
cooling  provided  by  electrically  powered  cryogenic  coolers.  Data  is  collected  by 
gamma  spectrometry  software  resident  on  dedicated  personal  computers  at 
each  station.  Calibration  of  the  detector  is  perfomed  using  a  Europium-152 
(i52Eu)  source  placed  a  various  angles  along  one  meter  radius  arc.  A  diagram 
of  a  station  is  shown  in  Figure  3. 


Fig  3:  Radiation  monitoring  station 


A  typical  spectrum  collected  during  the  study  is  shown  in  Figure  4.  The 
0.51 1  MeV  annihilation  radiation  peak  is  largely  due  to  positron  emitting 
radionuclides  released  from  LAMPF.  With  the  exception  of  ’^^Cs,  which  is  the 
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result  of  above  ground  weapons  testing,  all  ot  the  other  nuclides  identified  are 
naturally  occuring  (NCRP  1976).  The  code  READCHN  (APPENDIX  A)  was 
written  to  convert  the  binary  data  from  the  multichannel  analyzer  to  American 
Standards  Code  for  information  interchange  (ASCII)  format  for  analysis. 


ENERGY  (keV) 


Fig  4.  Typical  gamma  spectrum  obtained  during  monitoring  period 
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Stack  Monitoring 

The  radioactive  gases  created  in  the  targets  and  beam  stops  located  in 
experimental  areas  A,  B,  and  C  are  collected  and  piped  to  a  plenum  located 
at  the  base  of  the  FE-3  stack.  The  amount  of  activity  leaving  the  FE>3  stack  is 
measured  using  a  5  liter  Kanne  ionization  chamber  (Fig.  5)  (Engelke  and 
Israel  1973).  From  the  charge  collected  and  an  empirically  determined 
activity  per  charge  collected  ratio,  the  total  stack  gas  radionuclide  air 
concentrations  in  Becquerel  m-^  (Bq  m-^)  is  calculated.  The  release  rate  is 
then  calculated  by  multiplying  the  air  concentration  by  the  stack  flow.  Stack 
velocity  is  logged  by  an  anemometer  placed  halfway  up  the  stack.  The 
Isotopic  constituent  (Table  1 )  measurements  of  the  stack  gases  are  performed 
on  a  weekly  basis  using  a  germanium-lithium  detector  gamma  spectrometer 
and  a  pulse  height  unfolding  algorithm. 

Another  stack,  identified  as  the  FE-2,  is  located  on  the  site.  The  output  of 
this  stack  is  less  than  one  percent  of  the  FE-3  stack  and  was  disregarded  for 
purposes  of  this  study. 
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Fig.  5.  Radioactive  stack  gas  monitoring  system 


Table  1.  Isotopic  composition  of  exhaust  stack* 


Isotope 

Radioactive 

half  life 

Photons  emitted 
MeV  (intensity  (%)) 

Percent  of 
total  activity 
released 

16N 

7.12  s 

2.75(1) 

6.13(69) 

1.24 

IOC 

19.3  s 

0.511  (200) 

0.717  (100) 

3.90 

140 

1.2  m 

0.51 1  (200) 
2.313(99) 

1.57 

150 

2.0  m 

0.511  (200) 

62.29 

13N 

10.0  m 

0.51 1  (200) 

12.46 

11C 

20.4  m 

0.51 1  (200) 

18.23 

110m 

1.293  (99) 

0.31 

*  Personal  communication  (1992),  Simmonds,  S.  Los  Alamos  National 
Laboratory,  NM. 
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PROCEDURE 

Using  stack  effluent  activities  obtained  from  the  stack  monitoring  strip  chart 
recorders,  the  15  minute  interval  meteorological  data  from  the  meteorological 
tower  available  through  the  laboratory  common  file  system  (CFS),  and  the  data 
collected  on  an  hourly  basis  from  the  two  monitoring  stations  located  on  the 
north  side  of  the  Los  Alamos  Canyon,  the  measured  site  radionuclide 
concentrations  can  be  compared  to  the  predicted  radionuclide  air 
concentrations  from  the  CAP88-PC  code.  This  code  can  be  tested  using  the 
measured  meteorological  and  source  term  data  over  a  range  of  conditions. 

Meteorological  Data 

Since  the  meteorological  data  is  stored  in  a  format  which  is  not  amenable  to 
further  processing,  several  computer  codes  were  written  as  part  of  this  study  to 
aid  in  the  extraction  and  analysis  of  the  data.  APPENDIX  B  contains  the  listing 
for  the  program  METTRAN.  This  code  extracts  the  required  data  from  the 
weather  file  which  has  been  downloaded  from  the  central  laboratory  computer 
system  and  creates  a  new  file  with  the  extracted  data.  The  second  code, 
MET_AVG  (APPENDIX  C),  reads  the  data  file  created  by  METTRAN  and 
calculates  one  hour  averages  from  the  15  minute  data.  The  one  hour  data  can 
then  be  compared  with  the  source  term  data  for  that  interval  and  subsequent 
plume  travel  near  the  monitoring  stations. 

The  hourly  averaged  data  is  used  as  input  for  the  CAP88-PC  program. 
CAP88-PC  requires  that  meteorological  data  be  in  a  stability  array  format.  The 
stability  array  is  a  joint  frequency  distribution  of  wind  direction  and  wind  speed 
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for  each  stability  class.  The  stability  class  is  a  typing  scheme  developed  by 
Pasquill  (1^61)  using  the  letters  A  (extremely  unstable)  through  G  (moderately 
stable)  to  classify  the  turbulent  condition  of  the  atmosphere. 

Several  methods  may  be  used  to  derive  the  stability  class.  The  original 
work  by  Pasquill  (1961)  used  solar  insolation,  cloud  cover,  and  wind  speed  to 
determine  the  stability  class.  This  method  suffers  from  substantial 
uncertainties,  and  atmospheric  conditions  classified  as  stable  may  be  unstable 
and  vice  versa  (IAEA  1980).  The  temperature  lapse  rate  method  uses  the  bulk 
vertical  temperature  gradient  between  two  levels  in  the  atmosphere  to 
characterize  the  stability  class.  The  advantage  of  this  method  is  in  the 
simplicity  of  measurement.  While  it  characterizes  stable  conditions  reasonably 
well,  it  is  not  as  reliable  in  unstable  conditions  (IAEA  1980).  A  method  based 
on  both  temperature  lapse  rate  and  wind  speed  (Table  2)  is  more  sensitive 
since  it  includes  wind  speed  as  an  additional  variable  (Vogt  1977).  Additionally, 
the  conditions  under  which  this  typing  scheme  was  developed  (elevated  release 
and  rough  terrain)  more  closely  match  the  situation  being  modeled. 
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Calculation  of  Air  Concentration 

Each  station  was  setup  to  record  continuous  one  hour  spectra  between 
September  28  and  October  3, 1992.  Background  spectra  were  taken  during 
periods  when  the  accelerator  was  not  operating.  A  shield  was  placed  around 
the  lower  portion  of  the  detector  in  an  attempt  to  minimize  the  contribution  of 
radionuclides  present  in  the  ground  to  the  readings. 

Calibration  of  the  detector  was  accomplished  by  an  ^^^Eu  calibration 
source.  The  source  was  placed  at  one  meter  from  the  detector  and  spectra 
collected.  The  corrected  source  activity  was  3.94  xIO^  Bq  on  July  7,  1992. 
Table  3  lists  the  principle  i52Eu  photon  emissions  and  their  intensities. 

The  angular  dependence  of  the  detector  was  determined  using  the  ’“Eu 
calibration  source  at  15®  intervals  between  0®  and  90®  along  an  arc  placed  one 
meter  from  the  detector.  Figure  6  displays  the  response  of  the  detector  at 
various  angles  and  energies.  The  angular  dependence  is  most  severe  at  lower 
energies  and  at  60®,  but  is  reasonably  flat  at  energies  above  350  keV.  The 
response  as  a  function  of  energy  and  angle  was  determined  for  the  energy  of 
interest  (0.51 1  MeV)  by  interpolation. 


Table  3. 

i52Eu  data 

Energy 

Intensity 

(MeV) 

(%) 

0.122 

37 

0.245 

8 

0.344 

27 

0.799 

14 

0.965 

15 

1.087 

12 

1.113 

14 

1.408 

22 

16 


- • - 0.122  MeV  0.344  MeV  0.779  MeV 

_ _ _ 

Fig.  6.  Response  of  the  HpGe  detector  vs.  angle  and  energy 

Since  the  quantity  of  interest  is  the  concentration  of  radioactive  material  per 
unit  volume  (Bq  m-^),  the  calibration  of  the  detector  using  a  point  source  has  to 
be  related  to  the  irradiation  of  the  detector  from  a  passing  cloud.  Since  all  of 
the  nuclides  of  interest  are  positron  emitters  with  subsequent  annihilation 
photons,  the  dose  from  beta  radiation  is  neglected  in  the  subsequent 
discussion. 

For  clouds  with  dimensions  that  are  small  compared  to  the  range  of  the 
gamma  radiation,  a  calculation  of  the  gamma  dose  at  a  given  point  must  take 
into  account  the  radiation  received  from  various  parts  of  the  cloud.  The 
following  equation  is  suggested  by  Mealy  (Mealy  1968); 

0.1616ppgnEYQo(Qx/Qo)(Il+kl2) 

'  U 

where: 

yD  =  gamma  dose  (rads); 

\x  =  total  photon  absorption  coefficient  (m*’)  in  air; 
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Hen  =  energy  absorption  coefficient  for  photon  in  air  (m-’); 

Ey  =  average  photon  energy  emitted  at  each  disintegration 
(1.022  Mev  dis-1): 

Qx  =  amount  of  source  material  (Ci)  remaining  in  the  plume  after  a 
travel  distance  x; 

Qq  =  amount  of  source  material  (Ci)  released; 

U  =  average  wind  speed  (m/s); 

k  =  (P  -  4enV^len• 


and  Ij  are  defined  as; 


I.  = 


4(2;i)''2pa, 


n- 

00 


mr 


1  (m-rr  I 

,2 


2a 


-e 


[  (m+r)^  ] 

1  J 


dr*  dt 


(2) 


[(m- 


= 


iJ  "> 


2ai 


-e 


[  (tn+T)^  I 

1  J 


dr»dt 


(3) 


where: 

m  =  distance  from  the  center  of  the  clound  to  the  receptor  (m); 

0|  =  standard  distribution  of  the  material  in  the  cloud  (m); 

and  I2  do  not  have  exact  solutions.  They  are  evaluated  by  numerical 
integration  and  obtained  from  published  tables. 

If  the  dimensions  of  a  homogeneous  cloud  of  gamma-emitting  material  are 
large  compared  to  the  distance  that  the  gamma  rays  travel,  an  equilibrium 
condition  is  acheived.  Under  these  conditions  the  rate  of  energy  absorption  per 
unit  volume  is  equal  to  the  rate  of  energy  release  per  unit  volume.  For  a  semi- 
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infinite  uniform  containing  x  Bq  m-^,  the  gamma  dose  rate  (centiGray  s-i)  to 
tissue  at  the  cloud  center  is  given  by; 

YDi>=0.0093EYx(y)  (4) 

\vhere  pg/p  is  the  ratio  of  the  air  density  at  sea  level  to  air  density  at  the  point  of 
measurement.  For  Los  Alamos  the  ratio  is  1.3  (Bowen  1992). 

The  dose  to  tissue  can  also  be  calculated  by  the  relationships; 


Dy  =4^1 


Men 


V  p ; 

(O.SllMeV)(Cpp) 


where. 


(5) 

(6) 


T  =  energy  fluence  of  0.51 1  MeV  photons  crossing  the  detector 
(J  m-2) 

Met/P  =  •TI3SS  energy  attenuation  coefficient  for  0.51 1  MeV  photon  is  tissue 
(cm2/g) 

Cpp  =  number  of  counts  is  the  0.51 1  MeV  photopeak  of  the  detector 

Ej  =  average  intrinsic  photopeak  efficiency  of  the  detector  for  0.51 1  MeV 
photons 

The  average  photopeak  efficiency  can  be  determined  from  the  relationship 


Ei 


1  ^  ^TtCpp 

—  X  >  - — 

n  [‘=01 


(7) 


S 


where: 


The  total  number  of  photon  emitted  from  the  source  at  a  given 
energy 
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The  solid  angle  subtended  by  the  detector  at  6;  approximated  by  the 
frontal  area  of  the  detector  since  the  distance  between  the  source  is 
much  greater  (>103)  {^3^  the  frontal  area  of  the  detector 


n  = 


6  =  The  angle  between  the  calibration  source  and  vertical  axis  of  the 

detector 

Substituting  equations  (5)  and  (6)  into  equation  (4)  and  solving  for  x: 

X  =  ^3.69  X 10“^ Bq  x  m“^  x  cnt“*|Cpp]  (8) 

CAP88-PC  Models 

CAP88-PC  uses  a  modified  Gaussian  plume  equation  to  estimate  the 
average  dispersion.  The  Gaussian  plume  equation  is  the  mostly  common 
method  and  produces  acceptable  results  under  most  conditions. 


Plume  Rise 

Plume  rise  is  modeled  by  Rupp’s  (1948)  equation  for  momentum  dominated 
plume  rise.  The  plume  rise  is  added  to  the  actual  physical  stack  height,  h,  to 
determine  the  effective  stack  height,  H. 

Rupp's  equation  for  momentum  dominated  plumes  is; 


Ah  = 


1.5vd 

P 


(9) 


where: 


Ah  =  plume  rise  (m) 

V  =  effluent  stack  gas  velocity  (m  sec-’) 
d  =  inside  stack  diameter  (m) 
p  =  wind  velocity  (m  sec-’) 
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Plume  Dispersion 

Plume  dispersion  is  modeled  with  the  Gaussian  plume  equation  of  Pasquill 
(1961)  as  modified  by  Gifford  (1976): 


X  = 


27COyOxP 


gH/2(y/oy)2]^^l_1/2((z-H)/a2)2]  ^g[-1/2((z+H)/oz)2lj 


where: 


X  =  concentration  in  air  at  x  meters  downwind,  y  meters  crosswind, 
and  meters  above  ground  (Ci/m^) 

Q  =  release  rate  from  stack  (Ci  sec-’) 

li  =  wind  speed  (m  sec-’) 

Oy  =  horizontal  dispersion  coefficient  (m) 

Oj  =  vertical  dispersion  coefficient  (m) 

H  =  effective  stack  height  (m) 

y  =  crosswind  distance  (m) 

z  =  vertical  distance  (m) 

The  downwind  distance  x  comes  into  the  equation  through  Oy  and  a^,  which 
are  functions  of  x  as  well  as  the  Pasquill  atmospheric  stability  category 
applicable  during  emission  from  the  stack. 

The  equation  is  applied  to  ground-level  concentrations  in  air  at  the  plume 
centerline  by  setting  y  and  z  to  zero,  which  results  in: 


X  =  — (11) 

CAP88-PC  does  not  provide  output  of  the  concentration  at  a  specific  point; 
instead  the  ground-level  concentration  in  air  over  a  sector  of  22.5°  is 
approximated  by  the  expression: 
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Xave=fX  (12) 

where  f  is  the  integral  of  the  expression: 

e[-y2(y/ay)2]  (13) 

in  equation  (10)  from  a  value  of  y  equals  zero  to  infinity  divided  by  y^,  the  value 
of  y  at  the  edge  of  the  22.5°  sector,  which  is  the  value  of  the  downwind 
distance,  x,  multiplied  by  the  tangent  of  half  the  sector  angle.  The  expression 
is; 


/e 


dy 


ys 


(14) 


Solving  equation  (14)  and  substituting  into  equation  (10)  results  in  the 
equation  for  sector-averaged  ground  level  concentration  in  air; 


X  = 


Qg[-1/2(H/oz)2] 

0.1587l7txo2ji 


(15) 


Plume  Depletion 

The  only  mechanism  considered  for  plume  depletion  in  this  study  was 
radioactive  decay.  Depletion  is  accounted  for  by  substituting  a  reduced  release 
rate,  for  the  original  release  rate  Q  for  each  downwind  distance  (Slade 
1968).  The  ratio  of  the  reduced  release  rate  to  the  original  is  the  depletion 
fraction. 

The  depletion  fraction  for  radioactive  decay  is. 

Qi/Q  =  e-V  (16) 

where; 

Xr  =  effective  decay  constant  in  plume  (sec-i) 
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t  =  time  required  for  plume  travel  (sec) 

Dispersion  Coefficients 

Horizontal  and  vertical  dispersion  coefficients  (Oyand  a^)  used  for 
dispersion  calculations  in  CAP88-PC  and  for  depletion  fraction  determination 
are  taken  from  recommendations  by  G.A.  Briggs  of  the  Atmospheric  Turbulence 
and  Diffusion  Laboratory  at  Oak  Ridge,  Tennessee  (Moore  et  al.  1979).  The 
coefficients  are  different  functions  of  the  downwind  distance  x  for  each  Pasquill 
stability  category  for  open  country  conditions  (Table  4).  Figures  7  and  8  show 
the  plots  of  these  functions.  These  differ  significantly  from  power-function 
equation  suggested  by  Slade  (1968)  and  may  result  in  more  conservative 
(higher)  dose  estimates. 


Table  4.  Horizontal  and  vertical  dispersion  as  a  function  of  Pasquill  category 


Pasquill  Category 

Ov(m) 

o,(m) 

A 

0.22  x(1 +0.0001  x)-’^ 

0.20  X 

B 

0.16x(1+0.0001x)-’^ 

0.12  X 

C 

0.11  X  (1+0. 0001  x)-’4 

0.08x(1+0.0002x)-’^ 

D 

0.08  X  (1+0. 0001  x)-’^ 

0.06  x(1 +0.001 5x)-’^ 

E 

0.06  X  (1+0.0001  x)-’4 

0.03  X  (1+0.0003x)-i 

F 

0.04  X  (1+0.0001  x)-’4 

0.016  x(1+0.0003x)-i 

G 

calculated  by  subtracting  half  the  difference  between 
the  values  for  categories  E  and  F  from  the  value  for 
category  F 
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Downwind  Distance  (m) 

Fig.  7.  Horizontal  dispersion  as  a  function  of  downwind  distance  from  the 
source  and  Pasquill  stability  class 


Downwind  Distance  (m) 

Fig.  8.  Vertical  dispersion  as  a  function  of  downwind  distance  from  the  source 
and  Pasquill  stability  class  . 
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RESULTS  AND  DISCUSSION 


In  order  to  correlate  the  results  to  different  atmospheric  conditions,  they 
were  divided  into  atmospheric  stability  categories  A  and  B,  C  and  D,  and  E 
and  F.  The  data  were  also  divided  by  day  and  night  observations  and  wind 
speed  intervals. 

The  A  and  B  category  is  shown  in  Fig.  9.  The  linear  regression  equation 
for  this  data  is  y  =  1 .05 .  x ,  with  r^=  0.57.  Several  outlying  points  near  the 
origin  are  probably  due  to  varying  radiation  background  levels.  The  average 
difference  between  measured  and  predicted  values  is  -3.0  x  10-^  Bq  m-^.  In 
this  case  the  predicted  values  slightly  overestimate  the  measured  values. 


Measured  (Bq  rrr  3) 

Fig.  9.  Measured  vs.  predicted  values  for  stability  classes  A  and  B 


Categories  C  and  D  are  shown  in  Fig.  10.  The  data  shows  poor  correlated 
(r2=0.02).  The  average  difference  between  measured  and  predicted  values  is 
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3.0  X  10-3  Bq  m-3,  with  the  predicted  values  underestimating  the  measured 
values. 


Fig.  10.  Measured  vs.  predicted  values  for  stability  classes  C  and  D 


Stability  classes  E  and  F  are  shown  in  Fig.  1 1 .  The  model  generally  failed  to 
predict  most  of  the  measured  values  leading  to  the  conclusion  that  the  data  is 
uncorrelated  (r^=  0).  The  average  difference  between  measured  and  predicted 
values  in  this  case  is  2.4  x  10-^  Bq  m-3  This  result  is  probably  due  to  an 
underestimation  of  the  diffusive  power  of  the  atmosphere  and  terrain  by  the 
model  with  increasing  atmospheric  stability. 
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Fig.  11.  Measured  vs.  predicted  values  for  stability  classes  E  and  F 

The  data  is  divided  into  day  and  night  values  in  Figs.  12  and  13.  The  data  is 
correlates  for  daytime  measurements  (y  =  0.69,  r^  =  0.46).  The  average 
difference  between  measured  and  predicted  values  is  -  4.0  x  10-^  Bq  m-^.  For 
nighttime  measurements  no  correlation  was  found  (r^  =  0).  The  results  for  day 
and  night  follow  the  foregoing  discussion  for  stability  classes  since  unstable 
atmospheric  conditions  predominate  during  the  day  and  stable  (atmospheric 
inversions)  during  the  night. 
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Fig  12.  Daytime  measured  vs.  predicted  concentrations  (r2=0.46) 


0.07 

0.06 

'I  0.05 

g  0.04 

I  0.03 

1  0.02 
Q. 

0.01 
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■ 
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Best  fit  line 
y  =  -0.04  +.01 


0  0.02  0.04  0.06  0.08  0.1  0.12  0.14  0.16 

Measured  (Bq  m'^) 


Fig  13.  Nighttime  measured  vs.  predicted  concentrations  (r2=0.01). 


A  positive  correlation  between  wind  speed  and  the  measured 
concentrations  was  found,  r^  =  0.33  and  0.18,  for  the  center  and  west  stations, 
respectively.  Because  of  the  short  half-lives  of  the  radioisotopes  involved,  the 
higher  wind  speeds  will  result  in  a  higher  amount  of  activity  in  the  plume  when 
it  reaches  the  measurement  point.  This  leads  to  the  conclusion  that  higher 
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wind  speed  would  increase  the  exposure  at  the  boundary  although  they  tend 
to  also  cause  greater  dispersion  in  the  plume  (wind  speed  being  inversely 
proportional  to  concentration  in  the  Gaussian  plume  equation).  A  tabulation 
of  the  results  along  with  relevant  meteorological  data  is  shown  in  Table  5. 
The  windrose  for  the  survey  period  is  shown  in  Fig.  14. 
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CONCLUSIONS 

In  general,  the  model  slightly  overpredicts  (5%)  downwind  concentrations 
for  the  unstable  atmospheric  conditions  which  generally  occur  during  daytime. 
The  results  under  these  conditions  are  acceptable  considering  the 
nonuniformity  of  the  terrain,  variation  in  the  meteorological  conditions,  and  the 
uncertainties  associated  in  deriving  the  air  concentrations  from  the  gamma 
spectroscopy  data. 

For  nighttime  and  stable  atmospheric  conditions,  the  model  was  not  useful 
for  predicting  the  downwind  concentrations.  This  is  thought  to  be  due  to  an 
underestimation  of  the  dispersion  of  the  atmosphere  by  the  model  under 
these  conditions.  The  model  generally  did  not  predict  any  concentrations  in 
sectors  adjacent  to  the  sector  in  which  the  wind  was  primarily  blowing  for 
stable  atmospheric  conditions.  This  led  to  a  large  number  of  false  negatives 
in  the  north  sector  when  the  wind  was  primarily  blowing  from  the  south- 
southwest  (into  the  north-northwest  sector).  The  measurement  data  in  these 
cases  shows  measurable  elevations  above  background  in  the  north  sector. 

Nighttime  values  may  also  be  affected  by  the  lower  average  wind  speeds  at 
night.  Bowen  (1987)  found  winds  at  the  East  Gate  tower,  located  on  the  same 
side  of  the  canyon  as  the  radiation  monitoring  stations,  to  be  considerably 
higher  than  those  measured  at  the  LAMPF  tower.  If  the  true  wind  speed  is 
higher  than  that  measured  at  the  LAMPF  tower,  the  dovmwind  concentrations 
predicted  by  the  model  would  underestimate  the  measured  concentrations,  as 
is  the  case. 

it  is  difficult  to  draw  a  general  conclusion  as  to  whether  the  dose  to  an 
individual  would  be  over-  or  underestimated  from  the  foregoing  discussion.  It 
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depends  on  his  location,  the  predomininate  wind  patterns,  and  atmospheric 
stability.  If  data  is  collected  during  the  entirety  of  the  next  accelerator  run 
cycle,  the  greater  volume  of  data  will  allow  more  precision  in  the  comparison 
of  measured  and  predicted  values. 

Varying  gamma  radiation  background  levels  caused  significant 
uncertainties  in  the  data  for  certain  measurements.  It  would  be  advisable  to 
take  background  measurements  as  close  in  time  to  the  actual  measurements 
as  possible.  This  is  difficult  because  the  winds  blow  predominately  from  the 
south  and  southwest  and  there  is  no  control  over  the  source.  Statistical 
confidence  could  be  improved  with  longer  counting  times.  In  some  cases,  it 
was  difficult  to  define  a  peak  at  0.51 1  MeV,  even  when  the  wind  was  blowing 
into  the  sector  where  the  measurement  was  conducted. 
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FUTURE  WORK 

All  three  monitoring  stations  will  be  in  operation  during  the  next  run  cyde 
(May  -  October  1993).  Additionally,  weather  stations  will  be  installed  at  each 
monitoring  station,  improving  knowledge  of  the  micrometeorological 
conditions.  Planned  upgrades  to  the  stack  monitoring  system  include  logging 
the  output  activity  in  near  real  time.  Currently,  only  daily  output  is  logged. 
This  introduces  error  in  the  source  temi  calculations  for  periods  shorter  than 
one  day;  the  beam  current  and  hence  the  source  output  having  been  found  to 
vary  significantly  over  the  course  of  a  day.  These  upgrades  will  allow  greater 
data  confidence  and  should  be  in  place  for  the  entire  run  cycle  next  year. 

An  improved  method  of  deriving  the  air  concentration  from  the  gamma 
spectrum  data  could  be  developed.  If  data  on  average  plume  dimensions  for 
each  stability  class  were  known,  a  stability  class  dependent  conversion  factor 
could  be  developed.  A  laser  identification  and  ranging  system  (LIDAR)  may 
be  effective  in  providing  more  information  on  plume  dimensions.  Monte  Carlo 
techniques  could  then  be  applied  to  improve  the  accuracy  of  photon  fluence 
estimations  at  the  detectors.  Consideration  should  also  be  given  to  using  a 
more  powerful  dispersion  model  incorporating  local  terrain  features,  such  as 
the  Terrain  Responsive  Atmospheric  Code  (TRAC)  model  (Hogdin  1986). 
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APPENDIX  A 
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/********  MCS.C  **************************************************************>*****/ 
/********  C  PROGRAM  TO  CONVERT  MCS  DATA  TO  SINGLE  COLUMN  **/ 


/********  AND  REMOVE  HEADER  INFO  . . 

/*********15  JUL  92  . . >*«>..«**«*^«*.>*>»*.».**»»*.*>**>..*******/ 


#include  <stdio.h> 

#include  <string.h> 

#inciude  <stdlib.h> 

#include  <ctype.h> 

FILE  *fjn,  *f_out; 
int  data[5000]; 

main() 

{ 

char  infile[30],  outfile[30]; 
char  line[80],  mline[80]; 
int  c,d,  index,  dummychar; 
int  i,  nl=0,  nc=0,  channum; 

/*—  GET  INPUT  AND  OUTPUT  FILE  NAMES  ***********/ 
printfC'Program  MCS\nReads  multichannel  scaler  data,  strips 
header,and  reformats  data\n\n"); 

do 

{ 

printf("\n  Input  file  name?  "); 
gets(infile); 

} 

while  ((fjn  =  fopen(infile,"r"))  ==  NULL); 


do 

{ 
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printf("\n  Output  file  name?  "); 
gets(outfile); 

} 

while  ((f_out  =  fopen(outfile,"w"))  ==  NULL); 

index=0; 

do 

{ 

do 

{ 

/***  READ  FIRST  LINE  OF  FILE  AND  CHECK  FOR  END  OF  FILE  ***/ 

fgets(line,80,  fjn); 
if  (feof(f_in))  break; 

r  SKIP  OVER  HEADER  */ 

/**  FIRST  CHAR  OF  EVERY  DATA  LINE  CONTAINS  <ESC>  (027  IN  ACSIl), 
r  SO  IF  IT  DOESN’T  IT  MUST  BE  A  HEADER  LINE  **/ 

do 

{ 

fgets(line,80,fjn); 

}  while  (line[0]!=27); 
do 
{ 

line[0]=' ';  /**  write  over  <ESC>  character  **/ 

r  SCAN  LINE  AND  EXTRACT  DATA  ***/ 

sscanf(line,"%d  %d  %d  %d  %d  %d  %d  %d  %d", 

&channum,  &data[index],  Sidatalindex+I],  &data[index+2], 
&data[index+3],  &data[index+4],  &data[index+5], 
&data[index-i^],  &data(index+7]); 


if  (  channum  !=  index )  puts(line); 
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index  +=  8; 
fgets(line,80,fjn); 

}  while(line[0]==27): 

if  (ferror(f_in)) 

{ 

perror("Read  Error"); 
break; 

} 

if  (index  >  4088)  break; 

}while(!feof(fJn)); 

)while(!ferror(fJn)); 

/*  READY  TO  WRITE  OUT  DATA  */ 
r  INDEX  HOLDS  #  OF  DATA  POINTS  */ 

printfC'Number  of  data  points  is  %d\n", index); 


for  (i=0;  i<index;  i++) 
fprintf(f_out,  "%d\n",data[i]); 


} 


fclose(fjn); 

fclose(f_out); 
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APPENDIX  B 

A  LISTING  OF  THE  METTRAN  CODE 
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U - H 

II  METTRAN.CPP  // 

//  C++  PROGRAM  TO  EXTRACT  TA-53  METEROLOGICAL  DATA  // 
//  FROM  EM-8  FILES  // 

//  7  AUG  92  // 

//  MODIFIED  30  SEP  92  // 

n - n 


#include  <fstream.h> 

#inclucle  <stdlib.h> 

#include  <time.h> 

#include  <stdio.h> 

#include  <conio.h> 

//**  CHARACTER  VARIABLES  TO  HOLD  INPUT  STRINGS  **// 

char  id1[2],  jday[4],  time1[5],  dum1[4],  dir3(4],  dir2[4],  dir1[4],  dum2[4], 
sddjr3[4],  sddir2[4],  sddjr1[4],  dum3[4],  spd3[4],  spd2[4],  spd1[4], 
dum4[4],  sdspd3[4],  sdspd2[4],  sdspd1[4],  dunn5[4],  w3[4],  w2[4],  w1[4], 
dum6[4],  sdw3[4],  sdw2[4].  sdw1[4l .  (ine_2[100].  Iine_3[100],  line_4[100], 
temp3[4],  temp2[4],  tempi [4],  temp0[4],  dummy[80]; 

int  linejndex  =  1 ,  count  =  1 ; 

//  **  FLOATING  PNT  VARS  TO  HOLD  CONVERTED  CHAR  STRINGS  **// 

float  Jdl ,  Jday,  _time1 ,  _dum1 ,  _dir3,  _dir2,  _dir1 , 

_sddir3,  _sddir2,  _sddir1 ,  _spd3, 

_spd2,  _spd1,  _sdspd3,  _sdspd2,  _sdspd1, 

_w3,  _w2,  _w1 ,  _sdw3,  _sdw2,  _sdw1 , 

_temp3,  _temp2,  _temp1 ,  _temp0; 


int  main(int  argc,  char  *argv[]) 
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//**  PASS  INPUT.  OUTPUT  FILES  AS  PART  OF  COMMAND  LINE  **// 

{ 

if  (argc  ==  1) 

//**  QUIT  IF  NO  ARGUMENTS  **// 

{ 

cerr «  "ERROR  ==  >  Usage;  Mettran  {input  file}  (output  file}\n"; 
exit(1); 

} 

int  fjn  =  1 ,  f_out  =  2; 

ifstream  fin(argv[fjn],  iosi.noCTeate); 

//**  ERROR  IF  INPUT  FILE  DOESN’T  EXIST  **// 

if  (fin.failO) 

{ 

cerr « "ERROR  ==>  input  file  does  not  exist  \n"; 
exit(1); 

} 

ofstream  fout(argv[f_out],  ios::trunc /‘discard  contents  if  file  exists*/) 

//**  ERROR  WITH  OUTPUT  FILE  **// 

if  (fout.fail()) 

{ 

cerr  «"ERROR  ==>  output  file  errorVi"; 
exit(1); 

} 

//**  OUTPUT  FILE  HEADER  **// 

//**  APPEND  FILE  SPECS  AND  PROCESSING  TIME  TO  HEADER  **// 
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timejltime;  //**  TIME  FUNCS  DEFINED  IN  <TIME.H>  **// 
time(&ltime); 


fout «  "Input  FileAt"  «  argv[fjn]  «"\n" 
«  "Output  FileAt" «  argv[f_out]  «"\n" 
«  ctime(  &ltime  )  «  "\n"; 


fout  «  "jday"  «'At"  «  "time" 

«  "dir3"  «'At"  «  "dir2"  «"\t"  «"dir1" 

«'^t"  «  "sddir3"  «'At" «  "sddir2"  «'At" 

« "sddiM" « "\t"; 

fout  «  "spd3"  «"\t"  «  "spd2"  «‘At"  «  "spd1" 

«'^t"  «  "sdspd3"  «"\t"  «  "sdspd2"  «'At" 

«  "sdspdl"  «'At"  « '\v3"  «'At"  «  "w2"  «'^t" 

«  "w1"  «"\t"  «  "sdw3"  «'At"  « "sdw2"  «'M" 

« "sdw1"  «"\t": 

fout  «  "temp3"  «"\t"  «"temp2"  «'At" 

« "tempi"  «'V  « "tempo"  «’An"; 

//  **  READ  LINE  FROM  INPUT  FILE  **// 

//  **  NOTE:  BECAUSE  OF  FILE  FORMAT,  THE  DATA  MUST  BE  READ  **// 
ir  AS  A  CHAR  STRING  FIRST  **// 

II  **  NOTE:  GET()  APPENDS  TERMINATION  CHARACTER,  SO  MAKE  IT**// 
//**  ONE  CHAR  BIGGER  **// 


while( !( fin.eofO  ) ) 

{ 

fin.get(id1,2).getGday,4).get(time1,5).get(dum1,4).get(dir3,4).get(dir2,4) 

.get(dir1,4).get(dum2,4).get(sddir3,4).get(sddir2,4).get(sddir1,4) 

.get(dum3,4).get(spd3,4).get(spd2,4).get(spd1,4).get(dum4,4): 


46 


fin.get(sdspd3,4).get(sdspd2.4).get(sdspd1,4).get{dum5,4).get(w3,4) 

get(w2.4).get(w1,4).get(dum6,4).get(sdw3,4).get(sdw2,4).get(sdw1.4) 

.ignore(1,7n'); 

fin.get(dummy,15).get(temp3,5).get(temp2,5).get(temp1,5).get(temp0,5) 

.ignore(51,'\n'); 


//**  SKIP  THIRD  DATA  LINE  **// 

fin.ignore(81,'\n'); 

//**  CONVERT  STRINGS  TO  FLOATING  POINT  NUMBERS  **// 

ir  USING  ATOFO  FUNCTION  FROM  STDLIB.H  **// 

ir  AND  ADJUST  FOR  DECIMAL  PLACE  IF  NEEDED  **// 


Jd1  =  atof  (id1); 

Jday  =  atof  Gday);  //  JULIAN  DATE 

_time1  =  atof  (time1 );  //  TIME 


.dir3  =  atof  (dir3); 
dir2  =  atof  (dir2); 
dir1  =  atof  (dir1); 
sddir3  =  atof  (sddir3) 
sddir2  =  atof  (sddir2) 
sddirl  =  atof  (sddirl ); 


//  WIND  DIR  @  46  M  AGL 
//  WIND  DIR  @  23  M 
//WIND  DIR@  11.5M 
//  STD  DEV  OF  DIR 


,spd3=  0. 1  *  atof  (spd3);  //  WIND  SPEEDS 
spd2=  0.1  *  atof  (spd2); 
spd1=  0.1  ‘atof  (spdl): 

sdspd3  =  0.1*  atof  (sdspd3);  //  STD  DEV  OF  WIND  SPD 
sdspd2  =  0.1  *  atof  (sdspd2); 
sdspdl  =  0.1  *  atof  (sdspdl ); 


w3=  0.1  *atof(w3); 
.w2=  0.1  *atof  (w2); 
.w1  =0.1*  atof  (w1 ); 


//  VERTICAL  WIND  SPEEDS 
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_sclw3=  0.01  *  atof  (sclw3);  //  STD  DEV  VERT  WND  SPDS 
_sclw2=  0.01  *  atof  (sclw2); 

_sdw1  =  0.01  *  atof  (sdwl ); 

_temp3  =  0.1  *  atof  (temp3);  II  TEMPERATURES 

_temp2  =  0.1  *  atof  (temp2); 

_temp1  =0.1  *  atof  (tempi ); 

_temp0  =  0.1  *  atof  (tempO); 

//**  OUTPUT  DATA  TO  USER  SPECIFIED  FILE  **// 

fout  «  Jday  «'¥' «  _time1  «"\t" 

«  _clir3  «'¥' «  _dir2  «'^t" «  _dir1 
«■^t"  «  _sddir3  «'M"  «  _sddir2  «'¥’ 

«  _sddir1  «  "\t"; 

fout «  _spd3  «'V  «  _spd2  «'^t’' «  _spd1 
«'V  «  _sdspd3  «"\t" «  _sdspd2  «"\t" 

«  _sdspd1  «'M"«  _w3  «"\t"  «  _w2  «'V 
«  _w1  «'^t"  «  _sdw3  «'^t"  «  _sdw2  «'^t" 

«  _sdw1  «'V"  «  _temp3  «'^t"«  _temp2  «’^t" 

«  _temp1  «"\t"  «  JempO  «' V; 

//**  UPDATE  INDEX  VARIABLES  **// 

linejndex  =  linejndex  +  3; 
count  =  count  +  1 ; 

} 

ir  CLOSE  FILES  **// 

_fcloseall( ): 
return  0; 


} 
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APPENDIX  C 

A  LISTING  OF  MET_AVG  CODE 


II**  C++  PROGRAM  TO  AVERAGE  15  MIN  WEATHER  DATA  OVER  *// 


//**  ONE  HOUR  INTERVALS  *// 

//**  1  OCT  92  *// 

//**“******** . * . * . . . ^.^11 

#include  <fstream.h> 

#include  <stdlib.h> 

#include  <time.h> 

#include  <stdio.h> 


float  Jday[5],  _time1[5],  _dir3[5],  _dir2[5],  _dir1[5], 

_sddir3[5],  _sddir2[5].  _sddir1[5].  _spd3[5], 

_spd2I5],  _spd1I5],  _sdspd3[51,  _sdspd2[51. 

_sdspd1t5],  _w3[5],  _w2[51,  _w1[51,  _sdw3[5]. 

_sdw2[5],  _sdw1[5],  _temp3[5],  _temp2[5],  _temp1[5], 
_temp0[5]: 

//**  PASS  INPUT.  OUTPUT  FILES  AS  PART  OF  COMMANDI  LINE  **// 
int  main(int  argc,  char  *argv[]) 

{ 

if  (argc  ==  1 ) 

//**  QUIT  IF  NO  ARGUMENTS  **// 

{ 

cerr «  "ERROR  ==  >  Usage;  Met_avg  {input_file} 
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} 

int  fjn  =  1 ,  f_out  =  2; 

ifstream  fin(argv[fjn],  ios::nocreate); 

//**  ERROR  IF  INPUT  FILE  DOESN’T  EXIST**// 
if  (fin.failO) 

{ 

cerr «  "ERROR  ==>  input  file  does  not  exist  \n"; 
exit(1); 

} 

ofstream  fout(argv[f_out],  ios;;trunc  /*discard  contents  if  file  exists*/); 

//**  ERROR  IF  FILE  ALREADY  EXISTS  **// 

if  (fout.fail()) 

{ 

cerr  «"ERROR  ==>  output  file  already  exists"; 
exit(1); 

} 

//**  OUTPUT  FILE  HEADER  **// 

//**  APPEND  FILE  SPECS  AND  PROCESSING  TIME  TO  HEADER  **// 

time_t  Itime; 
time(&ltime); 

fout « "Input  File:\t" «  argv[fjn]  «"\n" 


«  "Output  File:\t"  «  argv[f_out]  «'^n" 

«  ctime(  &ltime ) « "\n\n"; 

tout « '^tTA-SS  Weather  Data  Averaged  Over  1  Hour  intervalsXn"; 

fout  « "jday"  «"\t" « "time"  «"\t" 

«  "dir3"  «’M"  «  "dir2"  «’\t"  «"dir1" 

«'M" «  "sddir3"  «'M" «  "sddir2"  «  "sddirl"  « '^t": 

fout «  "spd3"  «'M"  «  "spd2"  «’^t"  «  "spd1"  «"\t"  «  "sdspd3"  «' 
«  "sdspd2"  «’^t"  «  "sdspdl"  «'V"  « "w3"  «'¥'  «  "w2"  «'V" 
«  "wi  '•  «'^t"  « "sdw3"  «’M"  «  "sdw2"  «'^t"  « "sdw1 "  «'^t"; 

fout  « "temp3"  «'\\"  «"temp2"  «'M" 

«  "tempi"  «"\t"  «  "tempo"  «'^n"; 

ir  INITIALIZE  ARRAY  **// 

for  (int  j  =  0;  j  <  5;  j++) 

{ 

Jday[j]  =  0; 

_time1[j]  =  0; 

_dir3[j]  =  0; 

_dir2[jl  =  0; 

_dir1D]  =  0: 

_sddir3[j]  =  0; 

_sddir2[j]  =  0; 

_sddir1G]  =  0; 

_spd3[j]  =  0; 


_spd2[jl  =  0: 

_spd1ljl  =  0; 

_sdspd3[)]  =  0; 

_sdspd2[j]  =  0; 

_sdspd10]  =  0; 

_w3G]  =  0; 

_w201  =  0; 

_w1[j]  =  0; 

_sdw3[j]  =  0; 

_sdw2[j]  =  0; 

_sdw1[j]  =  0; 

Jemp3[j]  =  0; 

_temp2[j]  =  0; 

Jemp1  [j]  =  0; 

_temp0[j]  =  0; 

} 

//**  SKIP  OVER  6  HEADER  LINES 

for  (int  i=  0;  i  <  6;  i++)  fin.ignore(200,'\n'); 

while  ( !( fin.eofO  ) ) 

{ 


for  ( int  j  =  0;  j  <  4;  j++) 

{ 


fin  »  Jday[j]  »_time1[j] 
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»  _dir3[jl »  _clir2[j]  »  _clir1  [j] 

»_sddir30l  »_sddir2[j] 

»_sddir1[j] ; 

fin  »  _spd3[jl »  _spd2[j] »  _spd10] 

»  _sdspd30]  »  _sdspd2[j] 

»  _sdspd10]  »  _w3[jl  _w2[j] 

»  _w1[j]  »  _sdw3[j]  »  _sdw2[j]  »  _sdw1  [j]  » 
_temp3[j]»  _temp2[j]  »  _temp1[jl  »  _tempO[j] ; 

//**  STORE  RUNNING  SUM  IN  ARRAY  ELEMENT  4  **// 

_dir3[4]  =  _dir3[41  +  _dir30]: 

_dir2[4]  =  _dir2[4]  +  _dir20]: 

_dir1[4]  =  _dir1(4J +  _dir1[j]: 

_sddir3[4]  =  _sddir3[4]  +  _sddir3[j]; 

_sddir2[4]  =  _sddir2[4]  +  _sddir2[j]; 

_sddir1  [4]  =  _sddir1  [4]  +  _sddir1  [j] ; 

_spd3[4]  =  _spd3[4]  +  _spd3[j]; 

_spd2[4]  =  _spd2I4]  +  _spd2[j]; 

_spd1[4]  =  _spd1[4]  +_spd10]: 

_sdspd3[4]  =  _sdspd3[4]  +  _sdspd3[j]; 

_sdspd2[4]  =  _sdspd2[4]  +  _sdspd2[i]: 

_sdspd1  [4]  =  _sdspd1  [4]  +  _sdspd1  [j]; 

_w3t4]  =  _w3(4]+_w3[j]; 

_w2[4]  =  _w2[4]  +  _w2[j]: 

_w1  [4]  =  _w1  [4]  +  _w1  [j]; 
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_sclw3[4]  =  _sclw3[4]  +  _sdw301: 

_sdw2[4]  =  __sclw2[4]  +  _sclw2[j]; 

_sdw1  [4]  =  _sdw1  [4]  +  _sclw1 01; 

_temp3[41  =  _temp3[4]  +  Jemp30]: 

} 

fout «  Jday[3]  «'V  «  _time1[3] 

«"\t"  «  _clir3[4]/4  «"\t"  «  _dir2I4]/4  «  _dir1[4]/4  «"\t" 

«  _sddir3[4]/4  «'\r  «  _sddir2[4]/4  «'M"  «  _sddir1[4]/4  «'V"; 

fout  «  _spd3[4]/4  «’¥'  «  „spd2[4]/4  «'¥' «  _spd1[4]/4  «'V 
«  _sdspd3[4]/4  «'^t" «  _sdspd2t4]/4  «'M"  «  _sdspd1[4l/4  «'^t" 
«  _w3[4]/4; 

fout  «’¥' «_w2[4]/4  «’V  «  _w1[4l/4  «'¥’ «  _sdw3[4]/4  «’¥’ 
«  _sdw2[4]  «'¥’  «  _sdw1[4]/4  «"\t’’ «  _temp3[4]/4  «'¥’ 

«  _temp2[4]/4  «'¥' «  _temp1[4]/4  «'¥’ «  _temp0[4]/4  «  "\n" ; 

ir  RESET  ARRAY  ELEMENT  4  TO  0  FOR  NEXT  HOUR'S  DATA  *// 
j  =  4; 

JdayOl  =  0; 

_time10]  =  0; 

_dir30]  =  0; 

_dir20]  =  0; 

_dir10]  =  0; 

_sddir30]  =  0; 
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_sddir2[il  =  0; 
_sddir10]  =  0; 
_spd3[j]  =  0; 
_spd2[jl  =  0; 
_spd1Dl  =  0; 
_sdspd3D]  =  0; 
_sdspd20l  =  0; 
_sdspd1|j]  =  0; 
_w3D]  =  0; 
_w2D]  =  0; 
_w1D]  =  0; 
_sdw3[j]  =  0; 
_sdw2[j]  =  0; 
_sdw1  [j]  =  0; 
_temp3[j]  =  0; 
_temp2[j]  =  0; 
_temp1[j]  =  0; 
_temp0[j]  =  0; 

} 

fcloseallO; 


return  0; 

} 
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APPENDIX  D 

A  LISTING  OF  THE  READCHN  CODE 


57 


/********** . . . 

/**  READCHN.C:  PROGRAM  TO  READ  AND  CONVERT  **/ 
/**  HEADER  AND  CHANNEL  FROM  A  .CHN  FILE  **/ 

r  MODIFIED  FROM  MAESTRO  II  MANUAL  CODE  **/ 

/**17AUG92  **/ 

y*****>^***^»>>*>>.»»^>,>. . 


#include  <stdio.h> 

#include  <stdlib.h> 

#include  <process.h> 

#define  CHN  -1 

FILE  *fjn.  *f_out: 

r  PASS  INPUT,  OUT  FILES  AS  PART  OF  COMMAND  LINE  7 

int  main(int  argc,  char  *argv[]) 

{ 

char  accLtime[32];  /*  BUFFER  FOR  TIME  AND  DATE  7 
int  f_type; 

unsigned  int  chan_offset,  count,  mca_num,  num_chans,  num_writ,  segment; 
long  int  livetime,  realtime,  chan_data; 

if  (argc  ==  1  |  argc  ==  0)  /*  QUIT  IF  NO  ARGUMENTS  7 
{printf(  "ERROR  ==  >  Usage:  READCHN  {input_file} 
{outputfile}\n"); 

exit(1); 

} 

/**  ERROR  IF  INPUT  FILE  DOESN'T  EXIST  **/ 
if((fjn  =  fopen  (argv[1],  "rb"))==  NULL) 

{ 

printf  ( "ERROR  ==>  input  file  does  not  exist  \n"): 
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exit(1): 

} 

else  printfC'Input  file  %s  openedXn",  argv[1]); 

r  ERROR  IF  OUPUT  FILE  ALREADY  EXISTS  **/ 
if  (( f_out  =  fopen  (argv[2],  'W+‘’))==NULL) 

{ 

printfC'ERROR  ==>  output  file  already  exists\n"): 
exit(1): 

} 

else  printf("Output  file  %s  opened\n",  argv[2]); 

/**  READ  AND  OUTPUT  HEADER  INFORMATION  FROM  CHN  FILE  **/ 

fread(&f_type,sizeof(int),  1  ,fjn); 
if  (fjype  !=  CHN) 

{ 

printfC'Not  a  valid  file\n"); 
exit(1); 

} 

fread(&mca_num,sizeof(int),  1  ,fjn); 

fread(&segment,  sizeof(int),1  ,f_in); 

fread(acci_time+12,sizeof(char),2,fJn); 

fread(&realtime,sizeof(long).1,fJn); 

fread(&livetime,sizeof(long),1,fJn); 

fread(acq_time,sizeof(char),2,fJn); 

fread(acq_time+2,sizeof(char),3,fJn); 

fread(accLtime+5,sizeof(char),3,fJn); 

fread(acq_time+8,sizeof(char),2,fJn); 

fread(acq_time+10,sizeof(char),2,fjn); 

fread(&chan_offset,  sizeof(int),  1  ,fjn): 

fread(&num_chans,  sizeof(int),  1  ,fjn); 
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fprintf(f_out.  'TYPE  =  %4i  MCA  #  %2i  SEGMENT  # 
%3i\n",f_type,mca_num, segment); 
fprintf(f_out,  "realtime  =  %10li  SECONDS,  LIVETIME 
=  %10li  SECONDSVn",  realtime/50,  livetime/50); 

fprintf(f_out.  "DATA  COLLECTED  AT  "); 

fwrite(acq_time+8,sizeof(char),2,f_out); 

putc(':',f_out); 

fwrite(acq_time+10,sizeof(char),2,f_out); 

putc(':',f_out); 

fwrite(acq_time+12,sizeof(char),2,f_out); 

fprintf(f_out. "  ON  "); 

fwrite(acq_time,sizeof(char),2,f_out); 

putc('-',f_out); 

fwrite(acq_time+2,sizeof(char),3,f_out); 

putc('-',f_out); 

fwrite(acq_time+5,sizeof(char),2,f_out); 
fprintf(f_out,"\nSTARTING  CHANNEL  =  %6i.  NUMBER  OF 
CHANNELS  =  %6i\n\n",  chan_offset,num_chans); 

r  OUTPUT  CHANNEL  DATA  **/ 

fprintf(f_out, "Input  File;  %s\n\n",  argv[1]); 
fprintf(f_out, "CHANNEL  DATA\n"); 
for  (count  =  0;  count  <  num_chans;  cc  nt++) 

{ 

fprintf  ( f_out,  "%7i\t\t", count); 
fread  (&chan_data,sizeof(long),1  ,fjn); 
fprintf  ( f_out,  "%1 1l\n",chan_data); 

} 

fcloseallO; 


} 
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APPENDIX  E 

EXAMPLE  CAPP88-PC  OUTPUT 


61 


CAPP8  8-PC 
Version  1 .00 

Clean  Air  Act  Assessment  Package  - 1988 

WEATHER  DATA 
Non-Radon  Individual  Assessment 
Nov  22, 1992  4:58  pm 

Facility:  LAMPF 
Address: 

City:  LOS  ALAMOS 
State:  NM  Zip: 

Source  Category:  FE-3  Stack 
Source  Type:  Stack 
Emission  Year:  92 

Comments:  DOWNWIND  CONCENTRATIONS  FOR  JULIAN  DATE  276,  TIME 
2000  .2100hrs 

Dataset  Name:  276_21 

Dataset  Date:  Nov  22, 1992  4:58  pm 

Wind  File:  WNDFILES\276  21.WND 


62 


Nov  22. 1992  4:58  pm  WEATHER  Page  1 


HARMONIC  AVERAGE  WIND  SPEEDS  (WIND  TOWARDS) 


Pasquill  Stability  Class 

Wind 

Dir 

A 

B 

c 

D 

E 

F 

G 

Frequency 

N 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

NNW 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

NW 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

WNW 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

W 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

WSW 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

SW 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

SSW 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

S 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

SSE 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

SE 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

ESE 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

E 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

ENE 

0.000 

0.000 

0.000 

0.000 

6.055 

0.000 

0.000 

1.000 

NE 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

NNE 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 
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ARITHMETIC  AVERAGE  WIND  SPEEDS  (WIND  TOWARDS) 


Pasquill  Stability  Class 


Dir 

A 

B 

c 

D 

E 

F 

G 

N 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

NNW 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

NW 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

WNW 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

W 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

WSW 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

SW 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

SSW 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

S 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

SSE 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

SE 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

ESE 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

E 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

ENE 

0.000 

0.000 

0.000 

0.000 

6.302 

0.000 

0.000 

NE 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

NNE 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 
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Nov  22, 1992  4:58  pm  WEATHER  Page  2 

FREQUENCIES  OF  STABILITY  CLASSES  (WIND  TOWARDS) 


Pasquill  Stability  Class 


Dir 

A 

B 

c 

D 

E 

F 

G 

N 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

NNW 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

NW 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

WNW 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

W 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

WSW 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

SW 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

SSW 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

S 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

SSE 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

SE 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

ESE 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

E 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

ENE 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

NE 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

NNE 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

TOT 


0.0000  0.0000  0.0000  0.0000  1.0000  0.0000  0.0000 


ADDITIONAL  WEATHER  INFORMATION 


Average  Air  T emperature:  17.9  degrees  C 
291.1  K 

Precipitation:  0.0  cm/y 
Lid  Height:  1000  meters 
Surface  Roughness  Length:  0.010  meters 
Height  Of  Wind  Measurements:  1 0.0  meters 
Average  Wind  Speed;  6.302  m/s 


Vertical  Temperature  Gradients; 
STABILITY  E  0.073  k/m 
STABILITY  F  0.109  k/m 
STABILITY  G  0.146  k/m 


CAP88-PC  Version  1.00 
Clean  Air  Act  Assessment  Package1988 
CONCENTRATION  TABLES 
Non-Radon  Individual  Assessment  Nov  22, 1992  4:58  pm 


Facility:  LAMPF 
Address: 

City:  LOS  ALAMOS 

State:  NM  Zip: 


Source  Category:  FE-3  Stack  Source  Type:  Stack 

Emission  Year:  92 


Comments:  DOWNWIND  CONCENTRATIONS  FOR  JULIAN  DATE  276,  TIME  2000 
2100  hrs 


Dataset  Name:  276_21 

Dataset  Date:  Nov  22, 1992  4:58  pm 

Wind  File:  WNDFILES\276  21.WND 
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Nov  22.  1992  4:58  pm  CONCEN  Page  1 

ESTIMATED  RADIONUCLIDE  CONCENTRATIONS 
AT  VARIOUS  LOCATIONS  IN  THE  ENVIRONMENT 


Dry 

Wet 

Ground 

Air 

Deposition 

Deposition 

Deposition 

Wind 

Distance 

Concentration 

Rate 

Rate 

Rate 

Toward 

(meters)  Nuclide 

(pCi/m3) 

(pCi/m2/s) 

(pCi/m2/s) 

(pCi/m2/s) 

N 

700 

0-15 

O.OE+00 

O.OE-t-00 

O.OE+00 

O.OE+00 

N 

700 

N-13 

O.OE+00 

O.OE+OO 

O.OE+00 

O.OE+00 

N 

700 

C-11 

O.OE+00 

O.OE+00 

O.OE-t-00 

O.OE+00 

N 

743 

0-15 

O.OE+00 

O.OE-i-00 

O.OE+00 

O.OE+00 

N 

743 

N-13 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

N 

743 

C-11 

O.OE-i-00 

O.OE-^OO 

O.OE+00 

O.OE+00 

NNW 

700 

0-15 

O.OE+00 

O.OE-t-00 

O.OE+00 

O.OE+00 

NNW 

700 

N-13 

O.OE+00 

0.0E-»^00 

O.OE+00 

O.OE+00 

NNW 

700 

C-11 

O.OE+OO 

O.OE+00 

O.OE+00 

O.OE+00 

NNW 

743 

0-15 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

NNW 

743 

N-13 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

NNW 

743 

C-11 

O.OE-t-00 

O.OE+00 

O.OE+00 

O.OE+00 

NW 

700 

0-15 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

NW 

700 

N-13 

O.OE-i-00 

O.OE+OO 

O.OE+00 

O.OE+00 

NW 

700 

C-11 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

NW 

743 

0-15 

O.OE+00 

O.OE+OO 

O.OE-^00 

O.OE+00 

NW 

743 

N-13 

O.OE+00 

0.0E-»^00 

O.OE+00 

O.OE+00 

NW 

743 

C-11 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

WNW 

700 

0-15 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

WNW 

700 

N-13 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE-^00 

WNW 

700 

C-11 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

Dry 

Wet 

Ground 
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Wind 

Toward 

Distance 

(meters)  Nuclide 

Air 

Concentration 

(pCiym3) 

DeposKion 

Rate 

(pCi/m2/s) 

Deposition 

Rate 

(pCi/m2/s) 

Deposition 

Rate 

(pCi/m2/s) 

WNW 

743 

0-15 

O.OE+00 

O.OE+OO 

O.OE+00 

O.OE+00 

WNW 

743 

N-13 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

WNW 

743 

C-11 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

W 

700 

0-15 

O.OE+OO 

O.OE+00 

O.OE+00 

O.OE+00 

W 

700 

N-13 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

W 

700 

C-11 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

W 

743 

0-15 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

W 

743 

N-13 

O.OE+00 

O.OE-t-00 

O.OE+00 

O.OE+00 

W 

743 

C-11 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

WSW 

700 

0-15 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

WSW 

700 

N-13 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

WSW 

700 

C-11 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

WSW 

743 

0-15 

O.OE+00 

O.OE+00 

O.OE-i-00 

O.OE+00 

WSW 

743 

N-13 

O.OE-t-00 

O.OE+00 

O.OE+00 

O.OE+00 

WSW 

743 

C-11 

O.OE+00 

0  OE+00 

O.OE+00 

O.OE+00 

sw 

700 

0-15 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

sw 

700 

N-13 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

sw 

700 

C-11 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

sw 

743 

0-15 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

sw 

743 

N-13 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

sw 

743 

C-11 

O.OE+00 

O.OE-^00 

O.OE-^00 

O.OE+00 

ssw 

700 

0-15 

O.OE+00 

O.OE-i-00 

O.OE+OO 

O.OE+00 

ssw 

700 

N-13 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

ssw 

700 

C-11 

0.0E-»-00 

O.OE+00 

O.OE-^00 

O.OE+00 

ssw 

743 

0-15 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

ssw 

743 

N-13 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 
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Wind 

Toward 

Distance 

(meters)  Nuclide 

Air 

Concentration 

(pCi/m3) 

Dry 

Deposition 

Rate 

(pCi/m2/s) 

Wet 

DeposKion 

Rate 

(pCi/m2/s) 

Ground 

Deposition 

Rate 

(pCi/m2/s) 

SSW 

743 

C-11 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+OO 

S 

700 

0-15 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+OO 

S 

700 

N-13 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+OO 

S 

700 

C-11 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+OO 

S 

743 

0-15 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+OO 

S 

743 

N-13 

O.OE+00 

O.OE+00 

O.OE+OO 

O.OE+OO 

S 

743 

C-11 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+OO 

SSE 

700 

0-15 

O.OE+00 

O.OE+00 

O.OE+OO 

O.OE+OO 

SSE 

700 

N-13 

O.OE+00 

O.OE+00 

O.OE+OO 

O.OE+OO 

SSE 

700 

C-11 

O.OE+00 

O.OE+00 

O.OE+OO 

O.OE+OO 

SSE 

743 

0-15 

O.OE+00 

O.OE+00 

O.OE+OO 

O.OE+OO 

SSE 

743 

N-13 

O.OE+00 

O.OE+00 

O.OE+OO 

O.OE+OO 

SSE 

743 

C-11 

O.OE+00 

O.OE+00 

O.OE+OO 

O.OE+OO 

SE 

700 

0-15 

O.OE+00 

O.OE+00 

O.OE+OO 

O.OE+OO 

SE 

700 

N-13 

O.OE+00 

O.OE+00 

O.OE+OO 

O.OE+OO 

SE 

700 

C-11 

O.OE+00 

O.OE+00 

O.OE+OO 

O.OE+OO 

SE 

743 

0-15 

O.OE+00 

O.OE+00 

O.OE+OO 

O.OE+OO 

SE 

743 

N-13 

O.OE+00 

O.OE+00 

O.OE+OO 

O.OE+OO 

SE 

743 

C-11 

O.OE+00 

O.OE+00 

O.OE+OO 

O.OE+OO 

ESE 

700 

0-15 

O.OE+00 

O.OE+00 

O.OE+OO 

O.OE+OO 

ESE 

700 

N-13 

O.OE+00 

O.OE+00 

O.OE+OO 

O.OE+OO 

ESE 

700 

C-11 

O.OE+00 

O.OE+00 

O.OE+OO 

O.OE+OO 

ESE 

743 

0-15 

O.OE+00 

O.OE+00 

O.OE+OO 

O.OE+OO 

Wind 

Toward 

Distance 

(meters)  Nuclide 

Air 

Concentration 

(pCi/m3) 

Dry 

Deposition 

Rate 

(pCi/m2/s) 

Wet 

Deposition 

Rate 

(pCi/m2/s) 

Ground 

Deposition 

Rate 

(pCi/m2/s) 

ESE 

743 

N-13 

O.OE-^00 

O.OE+00 

O.OE+00 

O.OE+00 

ESE 

743 

C-11 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

E 

700 

0-15 

O.OE+00 

O.OE-i-OO 

O.OE+00 

O.OE+00 

E 

700 

N-13 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

E 

700 

C-11 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

E 

743 

0-15 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

E 

743 

N-13 

O.OE-^00 

O.OE+00 

O.OE+00 

O.OE+00 

E 

743 

C-11 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

ENE 

700 

0-15 

1.9E+00 

3.4E-07 

O.OE+00 

3.4E-07 

ENE 

700 

N-13 

6.3E-01 

1.  IE-07 

O.OE+00 

1.  IE-07 

ENE 

700 

C-11 

9.8E-01 

1.8E-07 

O.OE+00 

1.8E-07 

ENE 

743 

0-15 

1.7E+00 

3.1E-07 

O.OE+00 

3.  IE-07 

ENE 

743 

N-13 

5.8E-01 

1.0E-07 

O.OE+00 

1.0E-07 

ENE 

743 

C-11 

9.1E-01 

1.6E-07 

O.OE+00 

1.6E-07 

NE 

700 

0-15 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

NE 

700 

N-13 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

NE 

700 

C-11 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

NE 

743 

0-15 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

NE 

743 

N-13 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

NE 

743 

C-11 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

NNE 

700 

0-15 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

NNE 

700 

N-13 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

NNE 

700 

C-11 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

NNE 

743 

0-15 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 
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Dry 

Wet 

Ground 

Air 

Deposition 

Deposition 

Deposition 

Wind  Distance 

Concentration 

Rate 

Rate 

Rate 

Toward  (meters)  Nuclide 

(pCi/m3) 

(pCi/m2/s) 

(pCi/m2/s) 

(pCi/m2/s) 

IE 

743 

N-13 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 

IE 

743 

C-11 

O.OE+00 

O.OE+00 

O.OE+00 

O.OE+00 
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